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LARCH ARHA TON CHAMBERS FOR A HIGH ENERGY COSMIC RAY EXPERIMENT 


I . 

N77-25048 

INTRODUCTION 


Ion chambers are a prime choice for selection as a dE/dx detector 
for cosmic ray charge identification. They exhibit good resolution and 
linearity over a wide energy deposition range without interposing much 
mass in the particle beam, and have reasonably uniform response over 
large detector areas. 

The pulse ion chambers described in this report originally formed 
part of the High Energy Cosmic Ray Experiment for the high energy astro- 
nomical observatory (HEAD), (Orines et . al . , 1970), and also part of 
the Marshall Space Flight Center - UAH cosmic ray experiment. The objec- 
tive of the balloon flight program is to perform high resolution measurement 
of the iron group nuclei with the intention of determining individual 
abundances in the cosmic ray flux at the earth of particles with z = 22 
to 30. The experience being gained with these types of detector system 
is allowing a study of the feasibility of building much larger area 

di rectors to extend the charge resolution measurements to ’ igher z 
and energy. 

The balloon flight system which has now flown successfully a total of 
three times is shown in Figure 1. Of the three types of dE/dx detectors 
used in flu* experiment, pulse ion chambers, Cerenkov radiators and 
scintillators, the first has the best z resolution over a wide range of 
particle z and energy. The plastic scintillators define the particle 
acceptance cone of the telescope while the Cerenkov detector aids in 
rejection of low energy particle background and measurement of particle 


ve loc i ty . 


BALLOON BORNE COSMIC-RAY EXPERIMENT 



FIGURE 1. 
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This report includes work performed by the MSFC-UAH cosmic ray group 
as a whole. T. A. Parnell (MSFC) has been intimately involved in the 
project from the ou tser; G. A. Guenther and U. Pollvogt, formerly with 
the UAH, were involved in the early stages of design; .1. Derrickson (MSFC) 
performed some of the gas-studies; and T. A. Rygg , NAS Research Associate, 
computed the estimated integrated particle distributions for the balloon 


flight system. 


THEORETICAL APPROACH 


When ' protonlike ' particles pass through matter they lose energy 
by ionization and excitation of absorber atoms, Cerenkov radiative 
loss, bremsstrahlung, and nuclear interactions. The electromagnetic 
radiation losses are negligible compared to the ionization and 
excitation energy losses. Nuclear interactions are rendered very 
improbable by the low mass thickness of the detectors; a particle event 
being rejected if such an interaction occurs. 

For charged particles with rest mass much greater than that of the 
electron, the energy loss by ionization in an absorbing medium is given 
by the theoretical expression attributed to Bethe and Bloch (fc.g. see 
Livingston and Bethe, 1937): 


- dE /„ 2 
3 — = mtt r 
dx o 
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whore me 2 = rest energy of electron (=.51 MeV) 

r 2 = classical electron radius, c 2 /mc 2 
o 

I. = Avogadro's number 

p, Z, and A are the density, atomic number 
and atomic weight of the absorber 
7. ~ atomic number of incident particle 

B = v/c where v is the velocity of the particle 
I = average excitation potential per electron 
and may be taken as 11.5 x 10“^ MeV. 


In a thin absorber H and and dK /dx are assumed constant and the energy 
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loss given by the Bethe formula is proportional to the path length in the 

detector x p ■*. For ideal gases p/A is constant so that calculated 
A 

energy losses in different gas fillings for the same particle charge and 
energy are directly proportional to the ratio of Z's for the two gases. 

The specific energy loss is also seen from the Bethe formula to depend on 
the particle parameters and B. For balloon flights from Palestine, 

Texas the relativistic rise in dE/dx for particles with energies between 
the geomagnetic cutoff energy and 10 GeV or so is about 25-30%. In this 
energy internal the flux falls off by a factor of 30 approximately. 

For thin detectors the energy loss by a beam of monoenergetic particles 
is subject to rather wide fluctuations. These fluctuations, called straggling, 
are caused by statistical variations in the finite number of events which 
occur in the energy loss process. As a result, the energy loss distribu- 
tion for monoenergetic fast charged particles is quite wide and skewed to 
the high energy side since the number of large energy loss events is small. 
Theoretical treatment of energy loss straggling in thin absorbers was first 
made by Landau (1944). The problem was further analyzed by Symon (1949) 
and Vavilov (1957). 

Gas-filled ionization chambers and other thin detectors exhibit energy 
distributions considerably narrower than that predicted from Vavilov 
distributions. It must be emphasised that these detectors measure energy 
deposit within the detector volume rather than energy lost by the particle 
in the same volume. Watts (1972) has calculated energy deposit distribution 
functions for particles with several values of z up to 26 (iron) in gas-filled 
Ionization chambers. His calculations treated the transport of energy by 
A-rays (fast electrons produced by high energy transfer events between the 
relativistic particle and absorber electrons) from the point of primary 
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energy loss by the particle and along the A-ray track. Both energy trans- 
ported out of the detector volume and energy transported into the detector 
from prior energy loss events In the counter diaphragm were considered. The 
energy deposit distribution is consequently a very complex function of 
detector composition and geometry and does not lend itself to easy relative 
calculation say for a change of gas or counter thickness. 

An important result of this work was the quantified explanation of 
energy deposit distributions narrower than the Vavilov energy loss distribu- 
tions. The variation of distribution maximum goes as z^ in accordance 
with Bethe, while the variation of distribution width has an approximate 
z dependence. 

Since high energy A-rays produced above the detector volume deposit 
more energy in the counter than high energy A-rays produced in the volume, 
reductions in the mass above i he ionization chambers should result in a 
reduction of the width of the energy deposit distribution. 

Energy deposit distributions for particles of 10 GeV per nucleon and 
charges around z = 26 normally incident upon an ion chamber were calculated 
by Watts for our geometries and with argon and xenon gas fills. These 
curves have been replotted using cosmic ray abundances given by Jui I usson , 
Meyer and Mueller (1972) and arc shown in figures ?, 1 and U. 

Numbers for particles of each /. are given for a balloon flight of 10 hours 
duration at Palest ine, Texas. 

Because of Increased energy deposit in xenon compared with argon, the 
resolution of the iron peak is hc»i*i with xenon for the same thickness. ib. 
improvement is not simply in t lie square root of tin* ratio of their atomir 
numbers, since t h«* A-rays effects negate some nl the resolution increase. 


COUNTS/CHANNEL FOR 10 HR FLIGHT 
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A calculation was not made lor our specific geomet ry of 2 gaps 
of H.4 cm. of xonon each, hut an inexact estimate yields a resolution 
ol hette than 5 per cent FWHM for a monoonerget i c iron beam at several 
CeV per nucleon. 

Superimposed on the energy deposition frequency distribution are 
further uncertainties the most important of which are: 1) angular un- 
certainty of the particle track defined !>v the proportional counter 
hodoscope, 2) relativistic rise of energy deposit with particle energy 
J) charge collection efficiency within the ion chamber as a function of 
position, counter lifetime etc., and 4) amplifier noise. These effects, 
with the exception of 2), are discussed in the next section. 
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I A MORATORY DESfON STUDIES 

Mecha nica l Constru ctio n 

Each of the two ion chambers is comprised of 3 stainless steel mesh 
planes, the outer two electrically insulated from the middle, and all 
three from the housing. The square housings of internal dimension 67.4 cm 
are of aluminum (chosen for its low density) 0.125 in thick; one chamber 
being welded out of plate and the other machined from a single piece. 

One bousing received a chromate conversion coating and the other was 
anodized. Since surface finish is of great significance in the lifetime 
expectancy of gas-filled detectors, it was hoped to make a comparison. 
However, for both balloon flights we have operated the two housings bolted 
together with a common Xe-7% CH/ fill. Further details of the housing 
preparation and vacuum integrity testing are included in appendix A. 

Three photographs of the electrode structure and housing are shown in 
Figures 5, 6 and 7. Figure 5 shows an electrode assembly prior to install- 
ation in a housing. Figure 6 a closeup of the electrode frames in the 
vicinity of the feedthroughs, and Figure 7 an overall view. The dark spot 
in the center of the top cathode is an Am-241 source. 





cfovri nviuujWii 
riiifflonarwiiciii 





FIGURE 6 
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Each mesh plane is tightly stretched over square aluminum frames, 

62 cm outside dimension and 52 cm inside. Two frames with a total thick- 
ness of 2 cm hold each plane rigidly, some metal being removed to reduce 
weight. The mesh scieens are of 0.003 in (76 micrometer) stainless steel 
wire, 100 wires in - *. (A0 mm~l) with an optical transparency of 50%. 

Spacing is maintained between anode and cathodes, and between plane 
assembly and the housing, with ceramic and plexiglas insulators. Organic 
contamination is not considered a problem in ion chambers (unlike propor- 
tional counters). The diaphragms are 0.25 in thick aluminum honeycomb 
with a solid frame and 0.020 in aluminum faceplates. The diaphragms are 
sealed to the housing with Viton "0" rings. 

There are two insulated feedthroughs for the anode and cathodes and 
one 1.5 in port to which is attached a Viton-seated high vacuum valve which 
remains with the counter after gas fill. 

A stainless steel foil 0.002 in thickness is placed between the two ion 
chambers to render the outputs independent of each other. 

F h* 1 d On i fo rmity Che ck s 

The gap width of ion chamber I after final assembly for a balloon 
flight was checked. Measured deviation from the nominal A. 2 cm (1.65A in) 
nowhere exceeded * .010 in or * 0.6%. At the normal operating voltage of 
1000 V this produces an electric field variation of ♦ 6 V ; and since the 
measured charge collection efficiency has a slope of only 0.85% per 100 V 
at this operating voltage, variation of charge collection efficiency due 
to mechanical variation in gap width is seen to be negligible, (1 0.05%). 

Field uniformity at the edges of the ion chamber grids was investigated 


l>y f ie 1 d-p I ot t i ng. 


The electrode structure was painted at a scale of 2:1 op 
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conducting paper ns shown in Figure 8 and isopotential lines were con- 
structed. It may be seen that at the edges of the active area the frames 
tend to pinch the isopotentials (increasing the field) in the middle of the 
gap and weaken the field in the corners. The two effects tend to cancel 
in their impact on charge collection efficiency. Estimated field purturh- 
ations are about * 12% at l cm from the frame and * 6% 2 cm from the frame. 
The accuracy of the plotting technique is not better than * 2%. Possible 
reduct ion in charge pulse height for a particle track 2 cm from the edge 
of tlie act j ve area is estimated to he no more than 0.25% and about twice 
this at 1 cm. from the edge. 

Angu lar Uncertainty 

The angle at which a particle passes through the stack of detectors 
is defined by t he proportional counter hodoscopcs. These are four x-y 
pairs in two sets ol two pairs, tin* two sets being separated by 50 cm. 

Wire spacing between individual anodes is 0.5 cm. Angular resolution is 
- 1° which produces negligible track length uncertainty in the ion chambers 
at angles close to the Instrument axis (normal to the detector planes) hut 
increasing to about 17 . at maximum angle between track and axis. 

Ion Chamber Oas Studios 

Although there is a clear advantage in charge resolution of relativ- 
istic cosmic rays in using, xenon as t lie chamber gas fill instead of argon, 
a series of tests were made with both these gases in various mixtures with 
such ga* es as C (>2 ** n< l f’H/ 4 commonly u: ed to moderate electron temperatures. 
The cylindrical test ion chamber used in these studies is shown in Figure 4 
The chamber Is constructed ol stainless steel and has internal dimensions 
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large proportion of the CO 2 remains in the tubing between the ('.C >2 regulator 
and the ion chamber, and erroneous readings will be obtained if one proceeds 
with the experiment. To prevent this problem from arising, the gas-fill 
line was made of 0.25 in diameter stainless steel tubing to minimise its 
volume, and the pressure in the ion chamber was lowered about 10% before 
admitting CO 2 . After adjusting the partial pressure of CO 2 in the ion 
chamber to the new value, the total pressure in the chamber was returned 
to 760 torr with xenon. This had the effect of sweeping any C0 2 in the inlet 
manifold into the chamber ensuring that the chamber contained the calculated 
proportions of the two gases. 

The distribution of pulse heights for 5.5 MeV a's in the gas mixture 
was then obtained at anode voltages from 0-1200V. The channel number of the 
center of the peak of these distributions was plotted against applied voltage 
These plots, together with pulse risetime as a function of anode voltage are 
shown for mixtures of Ar-CC^, Xe-C0 2 , and Xe-CH^ between 1 and 10% of the 
quench gas; (see Figures 11-21). The data for the argon mixtures is shown 
for various values of amplifier shaping time, and pulse height clipping is 
evident where the shaping time is too short. 


Gas Contamination Effects on Charge Collection Efficiency 

Ion chambers are normally operated in the plateau region, which name 
refers to the plot of collected charge versus voltage applied. In this 
region of several hundred volts, essentially all the electronic charge is 
collected within several microseconds. The height, slope and onset voltage 
of this plateau are effected by gas contaminants, and the effects appear 
more pronounced with Xe then with Ar since the electronic drift velocity is 
lower in Xe. Charge losses may be either by columnar recombination, or by 
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electron capture: the latter reducing the drift velocity of the charge to 
the anode so that the fast electronics reiects it. Columnar recombination 
between electrons and positive ions has not been adequately treated. Jaffe's 
treatment is of recombination between negative and positive ions and requires 
that electrons be rapidly captured before they leave the track core. Electron 
capture and columnar recombination may be differentiated by careful experi- 
mentation. 

The migration of electrons in pulse ion chambers is strongly dependent 
on ihe composition of the gas. Not only does the drift velocity or col- 
lection time vary with the inert gas used but it is critically affected by 
small concentrations of electronegative gases. The importance of quantify- 
ing and then controlling or predicting these effects is of great importance 
lor gas filled detectors designed for long space flight, such as those which 
form part of some HEAO experiments. An example of the effects of contamin- 
ation on the performance of the cylindrical test ion chamber over an extend- 


ed period Is shown in Figure 22. Tills behavior is worse than we have obtained 


with the large area balloon chambers. Limited data for these is shown in 
Figure 21. In this case a peak pulse height versus voltage curve was not 
obtained, merely the peak pulse height al 1000V. Complete voltage 
character istics are given for two dates f> months apart in Figure 24 and 25 
for the same xenon-methane gas fill. The same plot is shown for a new gas 
fill into the same ion chambers in Figures 2h and 27. The chambers were 
pumped out and leak checked between fills without exposing the interior 
surfaces to the atmosphere at any time. Subsequent data after April 1974 
show a slower deterioration than for the September 1973 gas fill. Two 
effects are noticeable. Firstly, the level of I lie plateau falls with the 
charge collection efficiency. This is not a serious problem except at low 
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7 . since intermittent calibrations should always be performed. Rroaden- 
ing of the distribution must occur, however, due to poorer statistics. 
Secondly, the plateau threshold moves to the right so that the operating 
voltage is no longer on it. In this case field variations within the 
active volume will also cause charge resolution degradation. 

Pu 1 se K i set imes 

Risel lines were measured at the Tennelec preamplifier output and 
constitute the time taken for the pulse to rise between 10% and 90% of 
its full value. The measurements were made from photographs of the 
oscilloscope screen and in some cases, as noted on the figures, by direct 
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inspection of the screen. There is excellent quantitative agreement 
between measured ri set imes in all three argon-C02 mixtures and those 


calculated from the drift velocities given by English and Hanna (1953). 
Agreement between the two data sets for xenon-CX^ mixtures is only 
qualitative, which is probably due to the short (i-track length in xenon. 





Pr essur e Dependence of Pulse Heigh t 

An experiment was conducted with the cylindrical test ion chamber to 
measure the gas pressure dependence of charge collection ef r ieiencv. The 
gas used was Ar-10% CO 2 with gap width of 4.2 cm and a collimated Am-241 
source. Other operating conditions and the results are shown in Figure 28. 

It may be seen that the charge pulse corresponding to an americium 
u-particle is fairly constant around and above 1 standard atmosphere. The 
pulse height varies ♦ 1% with a variation in pressure of t 4%. (or t 30 torr) . 


Since the pressure in the ion chambers is normally adjusted to within 


torr, 


*• 


pressure variation is not considered a problem for a charge collection. It 
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should be emphasised that this experiment is rather specific for the 
set-up used: 4.2 cm of argon and 5.S MeV (t-part teles; however, it may bo 
deduced that there should be no difficulty in this respect with minimum 
ionizing particles. 

El ectro n i cs 

-14 

Detection of charge pulses of the order of 10 coulomb deposited 
in large area parallel plate ion chambers with capacitances of several 
hundred pf requires a very low noise charge-sensitive amplifier. A low 
power device to meet these requirements has been designed and used 
successfully on balloon flights. The diagram is shown in Figure 29. The 
detector (ion chamber anode) is ac-coupled to the FET input allowing the 
use of a large feedback resistor in l he charge loop. The FET must be 
specially selected for low noise. 

The design requirements of good signal/noise ratio at z=6 (carbon) 
are exceeded by the this amplifier-detector combination. Measured 
capacitance* between anode and caLhodo in both ion chambers was 260 pf 
approximately. Rms noise was 7mV measured at the postamplifier output, 
or about 6100 rms electrons, giving a signal to noise ratio of better than 
4:1 at z=3. 

A relat i vist ic iron nucleus <IeV per nucleon) deposits 

about 44 MeV of energy in a normal track through the xenon-filled balloon 
flight ion chamber. This corresponds to 2.01 x 10^ ion pairs, or 
1.21 x JO" 1 1 coulombs of charge ot cither sign. The amplifier (gain 12) 
produces a pulse of about 4 V for this charge which corresponds to about 
channel 400 of the 1024 channel AIM'.. 
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Ca lihr at ion 

Laboratory and in-flight calibration Is made with Am-241 sources of 
5.5 MeV 't-particles. These have a range of 2.1 cm in xenon at 1 atmosphere 
and normal temperatures. The foil-supported sources are mounted in the 
center of the cathode screens on the outside of the gap, one to each flight 
chamber. These are uncoil iraated in the present version. The source in the 
cylindrical test ion chamber is collimated producing a much narrower pulse- 
size distribution. An example of the effect of colliraation on this dis- 
tribution is given for some collimators. In each case the foil source is 
placed behind a hole or holes drilled through the 1/8 in. thick aluminum 
cathode disc. The collimating holes were: 

1) 1/H in. diameter hole, 1/8 in. long (45° spread) 

2) 1/16 in. diameter hole, 1/8 in. long (26.5° spread) 

J) 15 holes 1/32 in. diameter, 1/16 in. long over circle 1/4 in. 
diameter (26.5° spread) 

The variations of pulse height distribution is given for these three 
cases in terms of the KWHM; 1) 15.6%, 2) 6.4%, 3) 6.6%. The peak of the 
pulse distribution fell at channel 339 for case; 1) compares! to channel 
114 lor case; 2) and the tail on the higher charge side of the distribution 
extended further in case 1). 

This result indicates that the charge deposit from Am-241 in our 
configuration is not strictly comparable in terms of pulse height with a 
charge deposit from a relativistic particle. However, such an absolute 


calibration is not necessary. 
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SCINTILLATING PLASTIC FIBERS AS LIGHT PIPES 
FOR A COSMIC RAY HODOSCOPE: 


Feasibility calculations and measured atten- 
uation characteristics. 
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1 . Objectives of this Study 


Some cosmic ray experiments of the next decade will attempt to 
measure abundances in the high /. portion of the charge spectrum where 
fluxes may be few particles m -2 sr" 1 yr -1 . Such experiments require 
large areas to obtain adequate exposure and thus the development of 
charged particle hodoscopes of several square meters is needed. 

A candidate hodoscopc would use arrays of scintillator fibres, 
followed by an image intensifier and imaging system such as that proposed 
for the x-ray shadowgraph (Ref. 1). A literature search was performed 
to ascertain the experience of other workers with hodoscopes using this 
or similar principles. Calculations were performed to determine the 
feasibility of candidate systems and some laboratory experiments were 
performed to attempt to check these numbers. 

2. Plastic Scintillator Properties 


The properties of NK 102 fibres which make them interesting for our 
application are (a) their scintillating property, (b) their uniformity, 

(c) their availability in long lengths of various diameters, (d) their 
light guiding properties, (e) their low absorption coefficient for their 
own scintillation light and (f) the possibility of obtaining clad fibres. 

The plastic scintillating fibre we are interested in this article 
Is NE 102 by Nuclear Enterprises, Inc. with the following specification: 

Diameter 1.0 mm 

Index of Refraction 1.58 (not clad) 

Critical Angle of Total Reflection » 39.2652° 

Specific Gravity - 1.032 

Light Produced ■ 4230^ (2.9 ev) 
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Energy cl opor. J ti on of dJ fferent Z numb or particles and the total 
n Ult> h cr of jghof on r. gen cr ated in NE 102 fibres 



We performed calculations for three different charged particles: 

minimum ionizing 76 Ec and 6 C, and u particles from Am 241 decay. The energy 
dK 

deposition dx of each one in the fibre is 870 Mev /gm for minimum ionizing 

?c Fe; 46.8 ^ ev /gm for minimum ionizing c C; and the a energy is 

5 Mev approximately. Only a part of this energy can be converted to 
o 

photons (4230A, 2.9 ev) in the plastic scintillating fibres due to the 
saturation effect and conversion efficiency. The total number of photons 
(Np) can be calculated by the formula. 

E_ 

Np *= Ep x n x n sat 

where Np = Number of photons 

E » Energy deposition 

Ep *= Photon energy 

H ° Conversion efficiency (0.02 in our case) 
f| = Saturation coefficient (0.2 in our case) 

Sell 

The saturation coefficient is a function of dE/dx in the scintillator. 

If the energy deposit per unit volume exceeds some value, the emission centers 

no longer respond linearly. For minimum ionizing particles in plastic 

scintillators, n _ becomes less than one at about Z * 16. It is significantly 
sat 

less than one for stopping a's (~ 0.2). In these calculations n 8a t is taken 
as 0.2 for all three particles considered here, and consequently the 6 C 
case is treated rather conservatively. 
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For a particles in 1 mm thick fibre 

Np ■ 5,000,000 x 0.02 x 0.2 ■ 6,700 photons 
3 

For C particles in 1 ram thick fibre 

Np “ x 0.02 x 0.2 = 6,240 photons 

For Fe particle in 1 ram thick fibre 

N P — 1 — j 1 x 0.02 x 0.2 = 116,000 photons 

4 . Light Guide Properties of Fibres 


If the particle-induced photons are evenly distributed in 4 tt 
radians inside the fibre, then only a small fraction of these photons 
inside the total reflection cone can be transmitted by the fibre. This 
fraction can be calculated as follows: 



Solid angle of the total reflection cone 
where dA ■ rd0x27irxSin 0 

A ■ J~dA - J ^ 1 2tt r 2 Sin 0 d0 
■ 2 tt r 2 Cos 0 


dA 
— T 


0-0u » 90°-39. 2652° 

in 


0-0 


- 2tt r 2 x (0.3572) 


The fraction (f) 


f 



2n r 2 (0.3572) 
■— r 

4ti 


0.178 


That Is, if the reflectivity Is 1002 and the absorption coefficient is zero 
only 0.178 of the photons can be transmitted to one end. This result is 
independent of fibre diameter* 


4Q 


For a particle, there will be 6,700 x 1.178 * 1,200 photons 
which can be transmitted to PM tubes. 

For C particle, there will be 6,240 x 0.178 ** 1,110 photons 
which can be transmitted to PM tubes. 

For Fc particle, there will be 116,000 x 0.178 = 20,648 photons 
which can be transmitted to PM tubes. 

But the reflectivity and absorption coefficient become important if 
the guide length is long, (large number of reflections and long optical 
path) . 

The maximum number of reflections N p 
l R 

o L (total length) ^ L 
$ D R £ D tan 

A c The maximum optical length S 

. _ s _ 

S “ Co£Te c x N R = ' sin T ' ( Tnde P cnde,lt of D ' 

Table 1 shows the light transmission efficiencies of 1 ram diameter. fibres 
of various lengths. Table 1.1 shows the effect of absorption only: table 1.2 
of reflection loss only and l.e gives the total transmission efficiency using 
reflectivities obtained from Ref. 2. 

Table 2 and 3 convert the light produced by the 3 types of particle to 
photon numbers ar.d photoelcctron numbers, respectively, both quantities 
refering to the condition at the photo cathode. 


50 


TABLE I Total Light Transmission E ffic i en cy for 1 mm Diameter Fibres 

of Different Kinds 


1-1 Transmission Efficiency Due to Absorption 


Fibre Length 

100 cm 

50cm 

10cm 

Max Optical Length 
(For Ray at Total Internal 
Reflection angle) 

158cm 

79cm 

15. 


Absorption Coefficient 

0.015cnf 1 

0.018cm" 1 

0.015cm" 1 

0.018cm" 1 

0.015cm" 1 


Light Left 

9.32 

5.8% 

30% 

24% 

79% 

75% 


1-2 Transmission Efficiency Due to Reflection Loss 


Fibre Length 

100cm 

50cm 

lGcra 

Max Number of Reflection 

1224 

612 

122 

Reflectivity 

99.9% 

99.5% 

99.3% 

99.9% 

99.5% 

99.3% 

99.9% 

99.5% 

99.3% 

Light Left 

29.3% 

0.2% 

0.02% 

54% 

4.6 % 

1.3% 

88% 

54% 

42 % 

L_— — 


1-3 Total Transmission Efficiency [Eff ( abs.) x E ff (ref 1 . ) x 0.1 76] 


Fibre A* 

Reflectivity 

Absorption 
Coe If Jc 1 ont 

100 cm 

r-- — 

50 cm 

10 cm 

99.9% 

0.015cm" 1 

4.8 x 10" 3 

2.9 x 10" ? 

1.2 x 10" 1 

Fibre B 

99.5% 

0.015cm" 1 

3.3 x 10" s 

2.5 x 10" 3 

7.6 x 10" : 

Fibre C 

99.3% 

0.018cm- 1 

2.1 x 10- 6 

5.5 x 10-“ 

5.6 x 10- : 

J 


* Fibre A, very Rood fibre 

Fibre B, good nun-clad fibre 
Fibre C, regular fibre (Ref. 2) 















































TABLE 2 Photons reaching PM Tube for Different Sources and Different Fibres 


SOURCE 

Photon 

induced 

Photon can 
be transmitted 

100 cm 

50 cm 

10 cm 

Fibre A* 

Fibre B* 

Fibre C* 

Fibre A 

B 

C 

A 

B 

C 

a 

6,700 

1,200 

32 

0.22 

1.37 x 10” 2 

1.94 x 10 2 

16.3 

3.56 

8.1 x 10 2 

4.9 x 10 2 

z 

3.56 x 10 

C 

6,240 

1,110 

30 

0.20 

1.3 x 10" 2 

1.79 x 10 2 

15.3 

3.46 

7.7 x 10 2 

H 

9 

3.49 x 10 

Fe 

1.16 x 10 s 

20,648 

5.6 x 10 2 

3.8 

0.23 

3.2 x 10 3 

2.8 x 10 2 


1.4 x 10 4 

8.7 x 10 3 

6.4 10 s 


* See foot note on Table 1-3 


TABLE 3 Photoelectrons produced at photocathode PAGE IS POT 1 ? 


I 

100 cm 

50 cm 

! 

10 cm i 

SOURCE 

Fibre A B 

1 

C 

A 

B 

- 

C 

A 

B 

C 

i 

a 

3.2 j 0.22 x 10" 1 

1.37 x 10" 3 

1.94 x 10 1 

\ 

1.63 

0.36 

8.1 x 10 l 

4.9 x 10 l 

3.56 x 10 l 

C 

3.0 

0.2 x 10" 1 

1.3 x 10" 3 

1,79 x 10 1 

1.5 

0.3^ 

7.7 x 10 l 

4.7 x 10 1 

3.49 x 10 l 

Fe 

5.6 x 10 

0.4 

0.2 x 10- 1 

3.2 x 10 2 

2.8 x 10 1 

6^.4 

1.4 x 10 3 

8.7 x 10 2 


6.4 x 10 
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5. Quantum Efficiency and Background of PM Tube 

Table 3 shows the number of photoelectrons induced for different 
fibres and different sources at different distances. A typical photocathode 
quantum efficiency of 10 % was employed. 

A typical photocathode may emit 5, 000/sec, cm 2 thermionic electrons 
at room temperature. This is the main portion of the background. In our 
experiment we applied a 1*2 inch PM tube (« 11 cm ? in surface area) , which 
will emit about 55,000 electrons at room temperature. There would be .. 55 

output pulses/scc (0.1% of the 55,000 electrons) (Ref 3, A) with pulse 
heights corresponding to at least the average height for nine electrons 
starting at the photocathode. 

This means that for expected count rates, events could not be observed 
which produced only 9 photoelectrons (90 photons at photocathode). The 
crossed area in Table 3 showed difficulty of observation. But it should 
be quite possible to observe high Z particles over a 50 cm length with any 
one of the fibres on Table 3. 

6 . Experimental Results 

In the laboratory, we attempted to obtain Am ot spectra 
within 50 cm using NE 102 plastic scintillating fibre (Nuclear Enterprises, 

Inc.) in the arrangement shown in Figure 1. According to the manufacturer's 
report, the light loss of NE 102 is only about 0.009/cm corresponding to an 
attenuation length 2.5m. But unfortunately, we cannot observe any further 
than 2 inches (Figures 2, 3, A, 5, 6 ). The ends of the fibres were cut and 
polished and epoxied to microcope cover glasses. These in turn were coupled 
to the PM tube with coupling compound. 




Light-tight box 


hDistancej PM tube 


ibre — 


Sourcc(Ain 24l) 


Figure 1 The Experiment Set-up 



SCOPE 


Ch26 (7600) 


Numbers on curve are PHA 
channel numbers with counts/ 
channel in parenthesis. 


Ch 140 (149) 


Y 

Ch 46 (25) 


Ch 314 (5) 


Figure 2 Am 241 Source right on the PM tube. 
Accumulation time 450 Sec. 

Ch 23 (7740) 


Ch 107 (380) 


\ / 

Y 

*42 ( 86 ) 


QJ» 215 (13) 


Figure 3 Am 241 Source 0. 5" away from PM tube 
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Figure 2 shows Am 241 a-counts up to channel 314 on the PHA with a 

photoelectron noise spectrum cut off at about channel 46. If this latter 

corresponds to ~ 9 photoelectrons, a pulse in channel 314 corresponds to 

~ 63 photoelectrons. We calculated that roughly 120 photoelectrons should 

be produced under these experimental conditions; an estimate well within 

limits of error caused by uncertainty in the numerical values used. 

It is thought that surface damage such as scratch as, cracks, attached 

dirt, etc., perhaps produced by bad handling or bad end-cutting, was having a 

critical effect on the reflection efficiency, resulting in the loss of most 
photons within a gas centimeters. 

TABLE 4 Shows if the reflectivity dropped down to 95%, very few 
photoelectrons will be produced. 


SOURCE 

Reflectivity 

Light left after 5 cm 

Photons on PMT 

Photoelectrons 

a 

90% 

0.16% xO. 178 = 0.028% 

1.94 

0.192 

a 

95% 

4.3/. y 0.178 = 0.765% 

52 

5.2 

a 

96% 

8.2% x 0.178 = 1.46% 

99.4 

9.9 

a 

97% 

15.6% x 0.178 = 2.78% 

187 

18.7 


TABLE 4 


Fibure 7 is a plot of relative light intensity vs. distance. From the 
relation I *» I Q (where p is the reflectivity and N the number of reflection 
we estimated that the reflectivity (p) for 2 different fibres was 95-97%. From 
Table 4, we see that for an uparticle interacting with the fibre 5 cm from the 
photomultiplier tube, only 5-19 photoel cctrons are produced at the photocathode. 
This is of the same magnitude as the tail of the thermal noise spectrum. 
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7. Conclusions 

Plastic fibre scintillators 5x5x5 cm 3 have been used (Ref. 10) as 
beam finders for 700 MeV protons but no reports have been found of much 
longer hodoscopes. Other systems are described in Ref. 11 by Reynolds and 
coworkers (Princeton), Caldwell and coworkers (MIT), and an Imperial College, 
London, group. 

Glass fibres tend to have and maintain better surface smoothnes and 
Flanagan £t al . (Ref. 11) report detectability of Z=1 minimum ionizing 
particles at 20-25 cm. in 1 ram. diameter glass fibres. Experimental details 
were not given. They are not readily obtainable in long lengths ( > 20 cm.). 

Recent progress on optical fibre transmission line techniques (Ref. 5), 

-.'an reduce the surface damage produced by cutting to the minimum. Gloge, 

Smith, Bisbee £t al . at Bell Laboratories presented a method and their special 
tool to cut and connect fibres with little loss (Ref. 6, 7, 8, 9). 

Experimental results given here agreed well with calculated values of 
the photomultiplier-plastic scintillator output for Americium - 241 a-particles. 
The fibres used in the laboratory were found to have internal reflectivities 
(95-97%) lower than the range of those reported in the literature probably 
because we took no special precautions in handling them. It was demonstrated 
that, for fibres of the lengths of interest for a charged particle hodoscope, 
reflection losses are dominant. 

Tabel 3 shows that for aim. transmission length a reflectivity of 99.9% 
is required. This is difficult (but not impossible) with unclad fibres but 
not difficult with clad fibres. Thus, with good quality fibres, transmission 
distances of over 1 m. should be readily obtained. 

Further measurements should be made with clad plastic fibres and glass 


fibres. 
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N77-&5050 

EVALUATION OF CLAD SCINTILLATING LIGHT PIPES 
INTRODUCTION 


Arrays of fibers made of scintillating material have been 

used as position-sensitive detectors or hodoscopes for beam-finding 

at ion accelerators. Under this contract we have previously evaluated 

unclad fibers of NE 102 (Nuclear Enterprises, San Carlos) for use in 

a hodoscope for a heavy cosmic ray experiment. Experiments were made 
241 

with a's from an Am source incident upon one end of the fibre, the 

other end being viewed with a photomultiplier tube. The scintillation 

light could not be detected in any of the fibers tested beyond about 

5 cm. These previous results and calculations of the effective useful 

lengths for detection of relativistic heavy ions were given in a report 

and are reproduced in this Final Report. 

Light transmitted down a transparent fiber suffers a maximum 

of N * n reflections per cm, where r is the radius of the fiber 

2r tan 6 r 

c , 


i 


and n the index of refraction of the material (0 = arc sin — J The 

V c 

transmission of a length L of the fiber, assuming it is optically 
clear, is given by t * where a is the reflection efficiency. 

This expression shows that the transmission of such fibers depends 
very sensitively on the reflection coefficient, which in turn is a 
function of the surface condition of the fiber. This finish is depen- 
dent on the manufacturing process and how the material is handled 
afterwards, conditions which are very difficult to measure and control. 
To obviate this problem fibers may be coated or clad with another 
material of lesser refractive index which serves to protect the surface 
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of the fiber. The attendant light loss caused by the reduction of 
solid angle included in the cone of total reflectance is usually 
negligible when compared to the gains from increase in reflectivity 
at the fiber/cladding Interface. Transmission lengths of miles are 
now reported for clad glass light pipes. However, the change from an 
air/PVT interface (0 C ■ 39°) to a plexiglass/PVT interface (0 C ■ 71°) 
results in a drop in solid angle content of the cone of permitted 
light transmission fron about 18% to 3% for undirected light. This 
factor is independent of the fiber radius. 

Since no clad fibres of scintillating material were avail- 
able commercially, we have fabricated some and tested them using the 
same method as before. 

Fiber Manufacture 

Two scintillator materials were obtained from Nuclear Enter- 
prises, NE 102 and NE 150. The first is a polyvinyltoluene matrix 
and the second polystyrene. Each contained a proprietary and undis- 
closed mixture of scintillation chemicals. NE 102 has a claimed con- 
version efficiency of about 2% while NE 150 (no longer available) is 
less efficient. The materials were crushed from block form into 
particles which would pass a 0.25 inch sieve by Gold Plastics of Los 
Angeles. 

The granulated material was heated and extruded into a 
polymethylmethacrylate (plexiglass) sleeve which comprised the billet 
from which the fibers were drawn. Fiber manufacture was performed 
by International Fiber Optics of San Diego. 
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The following amounts of material were received: 





Nominal 





Fiber Core 

Diameter 

Length 




Material 

(inches) 

(ft.) 




NE 150 

0.04 

75 




NE 150 

0.06 

150 




NE 102 

0.04 

200 




NE 102 

0.06 

120 




NE 102 

0.08 

100 


Geometrical Characteristics of Clad Fibers 





Geometrical 

measurements are made on a 

fifteen-foot 

sample 

of 

each 

type of fiber. 

A wide variety of imperfections were noted 

but 

we are not able to 

discriminate between different types or 

magni 

tudes of 

imperfection in terms of quantitative effect on light 

trans 

mission. 

Typical diameters measured are given below. 



TYPE OF FIBER 

TYPICAL 

DIAMETERS (IN. 

) 

1 . 

NE 

150 .04" 

Major Axis .033 

Minor Axis 

.034 

2. 

NE 

150 .06" 

Circular .061 


3. 

NE 

102 .04* 

Major Axis .034 

Minor Axis 

.030 

4. 

NE 

102 .06" 

Major Axis .051 

Minor Axis 

.035 

5. 

NE 

102 .08" 

Major Axis .076 

Minor Axis 

.073 


Other Imperfections : 

1. Uniform except for knots up to 0.06 Inch diameter every few feet. 

2. Very circular but occasional large knots up to 0.11 Inch diameter. 
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3. Slight flattening but no major faults visible. 

4. Pitted surface. 

5. Not smoothly circular or elliptical like most of the fibers. 
Ridges present where fibers apparently stuck to each other after 
forming. 

Optical Characteristics 

The optical characteristics of the NE 150 fibers were so 
poor as to render them useless. For the plexiglass-clad NE 102 fibers 
we have: 


N 102 " 

1.58 


N pl. 

1.49 


0 

c 

70.6° 

(Critical angle for combination) 

^scint." 

4230A 

(Pt.ak of scintillation spectrum) 


Fibers were cut and polished and placed In a dark box with 

one end grease-coupled to a 2 Inch photomultiplier and the other close 

to a source of a particles (ab^^ ) . Lengths used were the longest 

for which these signals could be observed. 

The largest (0.080 inch) diameter fiber showed the best 

performance but even this was not very good. The maximum distance at 
241 

which an Am a-produced scintillation could be observed using the 
same experimental arrangement given previously was between 5 and 6 
Inches. Previous best transmitted distance with the unclad fibers 
was about 2 Inches though this was with a fiber of half the diameter 


(1mm). 
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A plot of logarlthum of the transmitted intensity (t) versus 
transmitted distance (figure 8) shows from the relation 


I - I p 
o 


NL 


(t - I/Io) 


a value of 92-95% for the reflectivity. 

This is no better than that found for the unclad fibers but 
an important uncertainty exists in the value of N. We have used the 
maximum value of reflections (i.e. those at the total internal reflec- 
tion angle, through the mean number of reflections suffered by a photon 
will be less). 

Conclusion 


Specially manufactured clad scintillating plastic fibers 
were tested in the laboratory and found to be of poor quality. Only a 
small improvement of a factor of 2 in transmission was achieved by 
cladding in this Instance. Much better fibers could doubtless be 
produced but this would require our spending much effort in fiber 
manufacturing technology. 
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Summary 

Thin fiat slabs of scintillator are a useful 
means of measuring angular location of gamma-ray fluxes 
of astronomical Interest. A statistical estimate of 
position error has been made of two scintillator 
systems suitable for ganma-ray burst location from a 
balloon or satellite platform. A single rotating 
scintillator with associated flux monitor Is compared 
with a pair of stationary orthogonal scintillators. 
Position error for a strong burst Is of the order of a 
few arcmln If systematic errors are Ignored. 


Introduction 


Slabs of scintillator may be used as directional 
gamma ray detectors for a number of studies in astro- 
nomy Including measurement of cosmological gamma ray 
anisotropy, location of discrete sources, measurement 
of atmospheric anisotropy, and location of gamna ray 
bursts . 

We examine here only the last of these, and 
attempt to estimate the position error of scintillator 
systems of reasonable dimensions used to locate these 
transient phenomena. 


Five years ago the gamma ray bursts were dis- 
covered by detectors on Vela spacecraft. Since then 
some 40 bursts have been detected by at least 9 space- 
craft In addition to several Vela satellites. The 
bursts are a unique phenomenon with no parallels In 
other energy regions or time scales of astronomy. 
Compared to other astronomical sources of gamma rays 
the bursts have large flux, a hard spectrum and 
duration of less than one second to some tens of 
seconds. They have not yet been Identified with any 
optical or x-ray object. Crude position measurements 
with uncertainties of 5 # -10* are available for less 
than a dozen bursts and so far indicate an apparently 
random distribution. 


The next major step In the understanding of 
these extremely energetic processes requires their 
position location to a degree of accuracy which allows 
Identification with an object previously observed In 
some other spectral region. Candidate Instruments 
and techniques which have been proposed for this pur- 
pose, include long baseline timing,* the x-ray 
shadowgraph, ^ an array of one dimensional Dlcke 
cameras,^ and an active antl-collimator.^ We examine 
here an alternative experimental arrangement using the 
directional sensitivity exhibited by a thin flat 
scintillator to a parallel beam of ganma rays. The 
general case of response of such shaped scintillators 
has been treated In detail by a number of authors. 


most recently Trombka et^. al_. (1975). We have 
selected here a scintillator of dimensions suited to 
the purpose of locating the direction of the gamma ray 
bursts from a small satellite. 

We have examined two configurations in detail. 
The first consists of a planar scintillator rotating 
rapidly about an axis in Its own plane. A single 
rotating scintillator would locate gamma ray bursts on 
an arc. Because of the rapid temporal variability of 
the burst flux, a second stationary scintillator 
is required to normalize the flux. Tlje second con- 
figuration considered is an orthogonal pair of 
stationary scintillators of similar dimensions. Both 
two-detector systems give a line position and require 
a third scintillator to obtain a point position for 
the burst origin. 

A thin flat detector rotating about an axis per- 
pendicular to the flux vector will respond to that 
flux according to a rectified cosine function, 

| cos 9j. If the position of the flux source is not 
known the function may be fitted to obtain the offset 
of the source from some reference direction. We have 
performed this fitting for some simulated bursts 
allowing for (a) counting statistics of the quantized 
flux (b) high local and diffuse background (c) Intense 
variation In flux with time over very short periods. 

Photon count rates in the detectors were esti- 
mated from published burst spectra^ and from natural 
background and used to estimate the location position 
resolution for various burst strengths. This res- 
lution was compared with that obtained by fixed ortho- 
gonal detectors of similar size. 


Rotating Scintillator 

For a '•hin flat scintillator of arbitrary shape, 
area A cm 2 and thickness d cm, the number of counts 
registered is given by 

N(t,0,t) • I j f A | cos 0| (l-exp(-a(E)d sec 0)} 

• I(E,t)dE dO dt...(l), 

where 0 is the angle between the flux direction and the 
detect r slab-normal, I Is the flux intensity In pho- 
tons cm -2 s - ’ keV“‘ and a is the abeorptlon coefficient 
for the detector material. The spin rate was chosen to 
be 2 revolutions s -1 with counts stored in 1024 bins/* 
radians. Dwell time per angular bln Is 0.244 ms. 

Counts per angular bln are summed throughout the burst, 
e.g. for a b'.rst length of 20s, the count rate in a 
particular angular segment is sampled 80 times. Each 
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sum of samples for each angular position is stored in 
one of 1024 bins. A similar set of 1024 sums is stor- 
ed in a corresponding register for the stntionary de- 
tector. Thus the i** bln in the register of one detec- 
ts covers exactly the same set of time perlo' a as the 
x 1 bin in the second register. 


The statistical error in 9 . measured over part 
or all of the burst is given by 


> 2 (0)- ( 2 - f ) [n 2 0 2 +N 2 o 2 +N 2 0 2 +N 2 o1 I 

\ | i I L 1 * 2 1 B * 2 R ‘ 2 Bi J ' 


The burst count SP summed over the burst period 
for each bin is calculated from equation (1). To each 
of these is added the mean background B for the same 
period, which has been previously determined. To 
simulate experimental reality, this set of numbers 
(SP+B) is randomized according to Poisson statistics 
and is called SP(R). The same process la performed 
for the stationary detector using a different set of 
random numbers to give ST(R). 


To remove the effect of temporal variation of 
burst flux, the normalized flux (AA) is calculated 
from 


AA 


SP(R) - B 
ST (R) - B 


AA still retains the form of |cos 0|, see figure 1, 
and is fitted to the trial function K|cos 9| where K 
is a constant of normalization (figure 2). The least 
squares method then gives the burst direction, 9. 

This was determined ten times using different random 
number sets and the standard deviation of 6 determined 
for several burst strengths. 


Calculations were performed for a 23 cm x 23 cm 
x 1 cm CsI(T£) crystal at a representative energy of 
100 keV. In practice a 1 cm Csl slab is rather effi- 
cient at this energy and the deviation frcxn the 
cosine function is not great. Edge effects were 
ignored . 

Several burst strengths were tried and the 
results of this least-squads f»t is shown in Table 1 
The strongest bursts have a total energy^of -4 x 10 -1 * 
erg cm -1 and a total of '400 photons cm - burst - in 
the energy Interval 40-200 keV are received at the 
Earth. The burst time profile of the April 27, 1972 
burst as recorded^ by the Apollo 16, shown in figure 
3, was used as a model and the count rates scaled up 
or down as tequlred. 


A total estimated background count of 1 photon 
cm -2 sec" 1 between 40 and 200 keV for the local 
(albedo plus charged-particle-produced) and diffuse 
components was used. This is assumed isotropic and 
Interacts with both front and back surfaces of the 
detector. This flux is typical for a low inclination, 
low altitude orbit. 


Stationary Scintillator Pair 

Consider the response of two stationary ortho- 
gonal thin flat scintillators to a y-ray flux from a 
direction perpendicular to the line of intersection 
of their planes. If 6 is the angle between the flux 
vector and the normal to one of the scintillators, 


tan 


■> 

N, 


where R, and R 2 are the actual count rates measured 
for each scintillator, and Nj, 2 and B, l2 are the 
respective numbers of burst and background photons 
respectively for each detector. 


The position error was calculated for the same 
burst parameters used for the rotating system and is 
shown in the table. Two cases were considered fpr 
each burst type; the burst flux vector in the plane 
of one detector, and the flux vector at 45° to both 
detectors. 


Conclusion 


The estimates of angular position errors made 
here indicate that the location uncertanties for 
gamma-ray bursts could be as low as a few minutes of 
arc for strong bursts and the sizes of scintillators 
considered if systematic errors in neasured rates can 
be kept below statistical errors. It is noted that 
statistical errors fall below 1Z for the "medium" 
bursts and these sizes of detector. With respect to 
systematic errors, the rotating scintillator will 
have an advantage over two stationary planar detec- 
tors. A direct comparison of two detector 

responses and efficiencies is required for the two 
fixed scintillators whereas for the rotating scintil- 
lator a comparison is >.nly required to correct for 
temporal fluctuations in the burst flux. 
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AN X-RAY SHADOWGRAPH TO LOCATE TRANSIENT 
HIGH-ENERGY CELESTIAL SOURCES 

G. J. Fishman and T. A. Parnell 
Space Sciences Laboratory 
NASA/Marshall Space Flight Center 
Huntsville, Alabama 

T. A. Rygg* and J. C. Gregory 
University of Alabama 
Huntsville, Alabama 

A new technique has been developed to locate 
strong, transient x-ray sources such as the recently 
discovered gamma-ray bursts. The instrument, 
termed a "shadowgraph, " locates sources by detect- 
ing the x-ray shadow cast by a large occulting mask 
pattern on an imaging detector. Angular resolutions 
of from 2 to 10 arc minutes are obtainable while 
essentially full sky coverage is maintained. The 
optimum energy range of operation is between 20 keV 
and 100 keV. 

The high-efficiency x-ray imaging detectors, 
which make it possible to locate bursts with intensi- 
ties down to ~10 photons/cm^ -sec, are capable of 
detecting single 20-keV photons with a spatial reso- 
lution of -0. 2 mm. The detectors consist of an 
x-ray to optical conversion phosphor, a multistage 
image intensifier, and a CCD image readout. 

I. INTRODUCTION 

The relatively new field of x-ray astronomy has led to the discovery of 
entirely new classes of celestial objects and has provided us with new insights 
into the final evolutionary stages of stars and the origin of high-energy pro- 
cesses in the galaxy. 

4 > 

Dr. Rygg died on March 14, 1975. 



Among the more interesting phenomena are the recently discovered 
gamma-ray bursts (Ref. 1). These infrequent, transient events produce intense 
fluxes of x-rays and gamma-rays relative to the background. Position data is- 
poor, although the sources appear to have an isotropic distribution. Although 
many theories for the bursts have been proposed (Ref. 2), it appears that a sat- 
isfactory explanation may have to await accurate measurement of their positions 
so that they may be identified with a particular type of optical or radio object. 
Such a measurement presents a unique experimental problem. Because of their 
apparently random occurrence in space and their brief duration — as short as 
0. 1 sec — it is necessary to use a wide-field detector. Since observable bursts 
occur infrequently and apparently randomly in time, the detector must be in 
continuous operation above the atmosphere for long periods. 

The instrument we have developed for locating discrete, brief bursts of 
x-rays and gamma-rays is termed a "shadowgraph" and combines the features 
of wide field, high efficiency, and good angular resolution. Image readout by a 
CCD confers an advantageous combination of simplicity, integrating capability, 
and energy discrimination of the x-ray spatial detector compared with other 
position -sensitive x-ray detectors. The x-ray system has thus far been tested 
only by photographic means rather than by CCD readout. Since CCD technology 
is expected to improve greatly before the first operational shadowgraph, this 
paper will stress the experimental techniques and capabilities of the shadow- 
graph and the specifications of CCDs optimized for this application. 

II. PRINCIPLE OF OPERATION 

The shadowgraph is composed of a large dome -like occulting shadow mask 
and a number of x-ray imaging detectors located inside to detect an x-ray shadow 
cast through any part of the shadow mask (see Figure 1). A typical mask is 
made from photoctched 0. 5 -mm -thick tungsten plate. This thickness would pro- 
vide adequate opacity for x-rays with energies up to 150 keV. The solid angle 
subtended by the mask at the detectors determines the field of view of the instru- 
ment, which could, conceivably, cover 4 tt steradians. 

Since it is assumed that the brief bursts arrive from a point source at a 
great distance (Ref. 1), the parallel incident x-ray beam would produce a sharp 
shadow of the coded mask on the face of an x-ray image detector. The x-ray 
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shadow projection angle is uniquely determined simply by finding the area on the 
mask which produced the shadow. The source direction is then derived from 
this shadow projection angle and the orientation of the spacecraft. 

The angular resolution of the shadowgraph is A0 = Ax/d, where Ax is the 
spatial uncertainty of the image position and d is the distance between the image 
plane and the appropriate shadow mask section. There are several components 
which enter into Ax; among the most important are: (1) x-ray source photon 

statistics and background events, (2) flexing and thermal distortion of the shadow 
mask and image detector image system, and (3) uncertainty in the location of 
individual x-ray photon events. It is estimated that for a weak burst, with -400 
x-ray photons over a 100-cm detector area, the value for Ax will be of the 
order -1 mm. For a small Explorer -type satellite, an average value for d is 
~70 cm, yielding a location uncertainty of A0 = 1 mm/70 cm = 5 arc minutes. 
Spacecraft motion during observation, and attitude uncertainty will increase the 
location error somewhat. 

The uncertainty due to photon statistics and background events, Axj, has 
been examined in detail, both analytically and through computer simulations. 
Consider N g photons incident on a shadow mask with a transmission factor r 
(ratio of open area to total area) and a characteristic pattern cell dimension £. 
Also assume that there are background events randomly spaced over the 
image detector area and that it is required to locate the pattern with a statistical 
uncertainty of n ff standard deviations. Then 


Ax 


1 


n 2 £(N + N, ) 

<r ' s b' 



( 1 ) 


Figure 2 shows the results of a simulation of x-ray images from various 
fluxes of signal photons from a point source incident on the shadow mask pat- 
tern, with r = 0. 44. A portion of the pattern is shown at the top of the figure. 
The amount of random background or noise photons was also varied, being equal 
to 1/2, 1/5, or 1/10 of the number of signal photons, as indicated in the figure. 
A computer program was used to find the optimum fit between the total mask 
pattern and the discrete photon image pattern by shifting the image in small 
steps over the larger mask pattern and finding the minimum photon coverage. 


Assuming that the image area is 10 cm X 10 cm, the optimized fit from the 
computer routine for the case of 200 signal photons and ICO noise photons was 
less than 0. 6 mm, while the value obtained from Equation (1) is -0-. 5 mm for a 
2 a error. Note that at this low signal level, the mask pattern is not visually 
recognizable. This demonstrates that relatively little information is required 
to ascertain the location of the source by this technique, since the image detected 
is a portion of a known pattern. Naturally, it would be impossible to reconstruct 
a meaningful image from only 200 photon locations. 

Since the locating accuracy for weak x-ray or gamma-ray bursts is depen- 
dent mainly upon photon statistics, it is desirable to obtain the highest ratio of 
signal photons to background photons in a particular energy range. Among the 
considerations that determine the optimum energy range are: (1) the spectrum 

expected from the source object, (2) the spectra of various components of the 
background such as the diffuse x-ray background and secondary x-ray and 
gamma-ray fluxes produced by charged -particle interactions; and (3) the effi- 
ciency of the shadowgraph mask and image detector combination. A rather 
detailed analysis of the combined spectral effects yields an optimum x-ray 
energy range of -20 keV to -100 keV in which to make observations of gamma- 
ray bursts. 

IE. IMAGING DETECTORS 

The imaging detector is the most critical component for the shadowgraph 
and other x-ray multiplex methods such as the scatter -hole or Dicke camera 
(Refs. 3, 4). For this application, the essential characteristics are: (1) good 
spatial resolution irrespective of angle of incidence of the x-ray photon, (2) high 
detection efficiency over the energy range of interest, (3) large sensitive area, 
and (4) capability for observing a large range of burst intensities. 

After studying several alternatives, the phosphor -image intensifier sys- 
tem (Figure 3) was chosen as the most suitable wide-field image detector for 
the energy region up to 100 keV. In this system, the x-ray photon is first con- 
verted into optical photons in a relatively thick (0. 5 to 1.0 mm) phosphor. 

Present phosphors under study include rare-earth phosphors and cesium iodide. 
The energy conversion efficiency of these phosphors is in the range from 5% tc 
15%, so that a 20-keV photon will result in “10- 5 photons in the range 4000 A to 
5000 A. The optical signal having a spot size of -1 mm^ is intensified and 
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demagnified by a multistage, electrostatically focused image intensifies The 
image is then transferred by a lens onto a CCD for image readout. 

Overall image reduction is from 80 mm X 80 mm at the x-ray phosphor to 
about 6 mm X 6 mm at the CCD. Linear spatial resolution at the x-ray phosphor 
of 0. 5 mm thus requires a 160 X 160 element CCD array of area about 6 mm 
square. Grey level encoding of 4 bits would be sufficient to provide coarse 
energy resolution for each detected x-ray. The energy information would be 
used in the data analysis to eliminate sources of background such as dark cur- 
rent, high-energy gamma-ray Compton events, and charged particles. It is 
anticipated that the x-ray image integration time will be of the order of 1 sec, 
which may require operating the CCDs at temperatures below 0* C. Such tem- 
peratures would also reduce the dark current of the image intensifier s. The 
total image data expected from a single x-ray burst is several megabits. Pres- 
ent CCD technology requires that this data be stored in a nonvolatile, onboard 
digital storage device (possibly a CCD memory) for later transmission to a 
ground station. 

The pictures in Figure 4 are output images of an experimental setup in 
our laboratory to study properties of intensifiers and x-ray converter phosphors. 

Two three-stage image intensifiers (Varo models 8585 and 8605) were coupled 
by an f/1.4 lens. A 0. 5-mm-thick rare-earth phosphor was placed directly on 
the fiber optic input of the first intensifier (25-mm diameter), and the output of 
the second intensifier was photographed on Tri-X film with a 35 -mm single -lens 
reflex camera. Figure 4 shows three photographs obtained with this system: 
two with weak radioactive sources present and one with bac^g ^und only. Indi- 
vidual x-ray photons can easily be observed from Co*^ (122 keV) and Am^^ 

(60 keV). All three photographs were obtained under the same conditions, out- 
lined in Figure 4. The centroid of the brighter x-ray images can be determined 
to within ~ 0. 2 mm. 

It appears that the imaging requirements for an x-ray shadowgraph in 
terms of array size, resolution, and dynamic range can be met by current CCD 
technology. For a space-borne experiment, a CCD image readout offers great 
advantages of simplicity, small size, low power, low weight, and image posi- 
tional stability compared to any conventional TV tubes. Areas of potential con- 
cern include long-term gain stability, image degradation during integration, and 
reliability. 
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Figure 1. Shadowgraph experiment components (As shown, an x-ray burst 

incident from the left would cast an x-ray shadow on the forward- 
facing image detector. ) 
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Figure 2. Simulated shadowgraph images 
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Figure 3. X-ray image detector schematic (The first-stage intensifier is an 
80 -mm/20 -mm zoom -type intensifier. Additional stages are 
required to produce a suitable intensity for the CCD. It may be 
possible to improve the overall efficiency by coupling the CCD with 
fiber optics or even by locating the CCD internally in the last inten 
sifier stage. ) 
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